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ABSTRACT 

This study confirms MLS110213:022733+130617 as a new eclipsing polar. We performed 
optical spectroscopic, polarimetric and photometric follow-up of this variable source identi¬ 
fied by the Catalina Real Time Transient Survey. Using the mid-eclipse times, we estimated 
an orbital period of 3.787 h, which is above the orbital period gap of the cataclysmic vari¬ 
ables. There are nine other known polars with longer orbital periods, and only two of them 
are eclipsing. We identified high and low-brightness states and high polarization modulated 
with the orbital period. The spectra are typical of polars, with strong high ionization emission 
lines and inverted Balmer decrement. The He II 4686 A line is as strong as H [ 3 . We modelled 
the photometric and polarimetric bright-state light curves using the CYCLOPS code. Our mod¬ 
elling suggests an extended emitting region on the WD surface, with a mean temperature of 
9 keV and B in the range 18 - 33, although the possibility that it could be a two-pole accretor 
cannot yet be ruled out. The WD mass estimated from the shock temperature is 0.67 Mq. 
The derived parameters are consistent with the eclipse profile. The distance was estimated as 
406±54 pc using the Period-Luminosity-Colours method. MLS110213 populates a rare sub¬ 
group of polars, near the upper limit of the period distribution, important to understand the 
evolution of mCVs. 
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1 INTRODUCTION 

AM Her, also known as polars, are a sub-class of magnetic Cata¬ 
clysmic Variables (mCVs): binary systems, composed of a highly 
magnetised (B ~ 7 — 120 MG) white dwarf (WD) primary star and 
a low-mass red dwarf secondary, in which mass is transferred from 
the secondary to the primary star by Roche lobe overflow. In con¬ 
trast to non-magnetic CVs, in which the transferred mass spirals in 
an accretion disc towards the WD, in AM Her stars the accretion 
occurs through a magnetic accretion column along the magnetic 
field lines and the rotation of both stars are synchronised with the 


* Based on observations made at the Observatorio do Pico dos Dias, Brazil, 
operated by the Laboratorio Nacional de Astroffsica. 
f Based on observations obtained at the Southern Astrophysical Research 
(SOAR) telescope, which is a joint project of the Ministerio da Ciencia, 
Tecnologia e Inova§ao (MCTI) da Republica Federativa do Brasil, the U.S. 
National Optical Astronomy Observatory (NOAO), the University of North 
Carolina at Chapel Hill (UNC) and Michigan State University (MSU). 
j E-mail:karleyne@gmail.com 


orbital period. Intermediate polars (IPs) constitute the other class 
of mCVs. They are not synchronised and may have an accretion 
disc that is internally truncated by the magnetic field of the WD. 
The optical and near-infrared radiation from polars is dominated 
by cyclotron emission and the X-ray radiation is dominated by 
bremsstrahlung, both from the accretion column near the WD, the 
post-shock region. Cyclotron radiation is highly anisotropic and po¬ 
larized. Therefore AM Her are the stellar sources with the largest 
fraction of polarized light in the sky, reaching 50 per cent in optical 
bands. See lCroppeil (11990h for a review on polars. 

According to the standard model, the evolution of CVs is 
based on angular momentum loss that decreases the separation be¬ 
tween the stars and, consequently, the orbital period. Thus, the sys¬ 
tems evolve from longer to shorter periods. There are two known 
mechanisms to caus e the angular momentum loss in CVs: the mag¬ 
netic braking (MB - fVerbunt & Zwaanlll981l) . which dominates in 
longer period systems and is driven by the outflowing magnetised 
stellar wind of the secondary; and the emission of gravitational 
waves (GW), which prevails in the shorter period systems. Most 
CVs have periods ranging from 1.2 to 6 h, but there is a deficiency 
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in the number of systems with periods between 2 and 3 h, known 
as the period gap. The period gap is understood as the transition 
from the MB to the GW mechanisms, considering that the sec¬ 
ondary is out of equilibrium in the MB phase due to the mass loss. 
When the MB halts, the secondary contracts to its equilibrium ra¬ 
dius and detaches from the Roche lobe, interrupting the mass trans¬ 
fer and becoming too faint to be detected. The mass transfer is re¬ 
established when the GW mechanism shrinks the Roche Lobe back 
to contact with the secondary, at orbital pe riods around 2 h (e.g., 
Howell, Nelson, & Rappaportj 120011 : [Webbink & Wickramasinghe 
20021 ) . 

The period gap is not so conspicuous in polars as it is 
in non-magnetic CVs. The presence of a strong magnetic field 
in the primary could make the MB mechanism less efficient, 
which is attributed to the trapping of the stellar wind in the 
combined magnetosphere of both _stars in synchronised sys¬ 
tems llWebbink & Wickramasinghe!I20021 . and references therein). 
Therefore, evolution of polars is slow and the accretion rates are 
low in comparison with others CVs. 

Another difference between polars and non-magnetic CVs is 
the scarcity of long-period AM Her systems. On the contrary, IPs 
are pre dominatly found above the period gap. IChanmugam & Ravi 
d 1984h proposed that IPs evolve to polars as the orbital period de¬ 
creases as result of standard CV evolution. A problem with this 
scenario is the different distributions of magnetic field intensities 
in polars and IPs. ICummingl d2002l) proposed that the WD mag¬ 
netic field can be reduced by accretion, which could redeem that 
scenario. However. lZhang et al . ! ( 20091 ) performed new calculations 
considering non-spherical accretion and concluded that the mag¬ 
netic field is not easily modified by accretion. Presently itjs not yet 
clear what is the evolutionary relation betw een mCVs dNorton et al.l 
l2008l : lPretorius. Knigge. & Schwopel2013l . and references therein). 

The Catalina Real Time Transient Survey (CRTS: lDrake et al] 
120090 improved the number of discoveries of faint transient ob¬ 
jects. CRTS is composed by three separate surveys - the Catalina 
Schmidt Survey (CSS), the Mount Lemmon Survey (MLS) and the 
Siding Spring Survey (SSS) - that repeatedly scan a combined area 
of 30 000 deg 2 of the sky in both hemispheres, up to a limit of 
19-21 mag in V band. This survey has already found more than 
one thousand CV candidates, which are made public on their web- 
sittE MLSl 10213:022733+130617 (hereafter MLS110213) is one 
of those candidates: it is a Galactic source that is highly variable, 
ranging from 16.5 mag to 19.5 mag in the CRTS light curve in 
timescales of days. We selected it, as a target for detailed observa¬ 
tional follow-up, from a larger program of spectroscopic classifica¬ 
tion of mCV candidates (Oliveira et al., in preparation). 

As of June 201 4/1 114 polars have been identified 
l l Ritter & Kolbll2003Ll201 lT hereafter RK catalogue). Among them, 
30 systems have periods above the period gap and only 5 of 
them are eclipsing, being MLS 110213 one of those systems. 
MLSl 10213 is included in the RK catalogue based on reports in 
VSNElE. There are only two other known eclip sing p olars having 


longer periods, namely V895 Cen ( 

Howell et al. 1997) and V1309 

Ori (Staude, Schwope, & Schwarz 

200]]). Here we report the re- 


suits of the photometric, spectroscopic and polarimetric follow-up 
of MLS 110213 and the first modelling of the light and polarization 


curves. In Section 2, we present the data. In Section 3, we show 
the analysis of the light curve, of the polarimetric data and of the 
spectra. The modelling of the light and polarization curves using 
the CYCLOPS code is presented in Section 4. The eclipse fitting is 
explained in Section 5. In Section 6 we estimate the distance to 
MLSl 10213. We discuss our results and present our conclusi ons i n 
Section 7. Preliminary results were presented in lSilva et al.l J20 14h . 


2 OBSERVATIONS AND DATA REDUCTION 

Photometric and polarimetric observations of MLS 110213 were 
performed at the Observatorio Pico dos Dias (OPD - LNA/MCTI), 
located in southeast Brazil, using the 1.6-m Perkin-Elmer tele¬ 
scope, while the spectroscopic data were obtained on the SOAR 
Telescope, on Cerro Pachon, Chile. The journal of observations is 
presented in Table Q] The data acquisition and reduction are de¬ 
scribed in the following sections. 

2.1 Photometry 

Photometric data in Rc, Ic and white light (clear) bands were ob¬ 
tained employing a thin, back-illuminated 2048 x 2048 E2V CCD 
42-40 mounted in an IkonL camera. The exposure times were 
scaled according to the atmospheric conditions, and acurate tim¬ 
ings were provided by a GPS receiver. Bias and dome flat-field ex¬ 
posures were used for correction of the detector read-out noise and 
sensitivity using IRAI0 standard routines. The differential aperture 
photometry was performed with the DAOPHOT II package, adopt¬ 
ing as reference the star USNO B 1.0 31—0032240 (R = 15.38 mag; 
I = 15.23 mag). The conversion between R and Land olt’s Rc fo r 
this object indicates a difference of only 0.03 mag l iKidgeill2003l) . 
The resulting light curves in Rc, Ic and white light are shown in 
Fig.0 

2.2 Polarimetry 

We obtained Rc polarimetric data using a CCD cam- 
era modifie d by the polarimetric module, as described in 
ll Magalhaes et alJ Il996b . The data reduction followed the 
standard procedures using IRAF. The p olarizat ion w as c alcu¬ 
lated acco rding to iMaga lhaes. Benedetti & Roland! J19841) and 
iRodriguesTc ieslinski & Steined (1998) using the package PCCD- 
PAC Kjjperevrall2000l) and a set of IRAF routines developed by our 
groupj. Each set of 8 images is used to produce one measurement 
of linear and circular polarizations. We grouped the images in 
the following way: 1-8, 2-9, 3-10, and so on. Hence the polariza¬ 
tion points are not independent measurements. The polarization 
position angle correction to the equatorial reference system was 
performed using standard polarized stars. The polarization of 
the non-polarized standards is consistent with zero, therefore no 
instrumental polarization correction was applied. We could not 
determine the correct sign of the circular polarization - if positive 
or negative. But we could detect changes in this signal. The linear 
polarization is always a positive value, which introduces a bias 


4 IRAF is distributed by the National Optical Astronomy Observatories, 

1 http://crts.caltech.edu which are operated by the Association of Universities for Research in 

2 update RKcat7.22, 2014. at http://physics . open. ac. uk/RKcat/ Astronomy, Inc., under cooperative agreement with the National Science 

3 VSNET, at http ://www. kusastro . kyoto-u . ac . jp/vsnet/ Foundation. 

The vsnet-alerts 12847, 12855, 12884, 12888 are related to MLSl 10213. 5 http://www.das.inpe.br/~claudia.rodrigues/polarimetria/reducao. 
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Table 1. Log of observations of MLS 110213. 


Date 

Telescope 

Instrument 

Filter 

Exp. time (s) 

Duration (h) 

2011 August 14 

1.6 m P-E 

Caml+IkonL 

Ic 

60 

2.7 

2011 September 22 

1.6 m P-E 

Polarimeter + IkonL + A/4 

Rc 

90 

5.0 

2011 October 27 

1.6 m P-E 

Caml+IkonL 

Rc 

100 

5.1 

2011 October 28 

1.6 m P-E 

Caml+IkonL 

Rc 

60 

2.5 

2012 August 13 

1.6 m P-E 

Caml+IkonL 

Ic 

60 

3.3 

2012 September 10 

4.1 m SOAR 

Goodman HTS 

- 

1200 

1.0 

2013 August 23 

1.6 m P-E 

Caml+IkonL 

Clear 

10 

1.5 



Figure 1. OPD light curves of MLS 110213 folded with the ephemeris 
estimated in Section [TT] From top to bottom: Rc band, Ic band and white 
light. White light magnitudes are not calibrated. 


into any measurement: the measur ed polarization value is gre ater 
than the true polarization value dSimmons & Stewart! 11985h . In 
particular, measurements having P/up < 1.4 provides only 
upper limits to the real polarization value. The polarization of 
MLS110213 was bias-corrected following IVaillancourtl l2006t ). 
Figure [2] shows the resulting polarization curves for MLS110213, 
binned in 40 orbital phases. The lack of points around phase 0 
corresponds to the eclipse in which we had no signal to measure 
the polarization. 

The ordinary and extraordinary counts of the polarimetric data 
were summed to obtain the total counts of each object. Hence the 
polarimetric data allowed us to perform differential photometry and 
to obtain an additional MLS 110213 light curve, which is also plot¬ 
ted in Fig.Q] 


2.3 Spectroscopy 

Spectra of MLS 110213 were obtained in queue mode on Septem¬ 
ber 10, 2012 using the SOAR Telescope, Chile , with the Good¬ 
man High Throughput Spectrograph dClemens. Crain. & Andersonl 
120041 ) and a Fairchild 4096x4096 CCD with 15 micron/pixel. The 



0.0 0.5 1.0 1.5 2.0 

Orbital Phase 


Figure 2. Polarimetry of MLS 1 10213 in Rc band. The data points were 
grouped in 40 orbital phase bins. The lack of points around phase 0 corre¬ 
sponds to the eclipse in which we had no signal to measure the polarization. 
From top to bottom, the panels show the circular polarization, linear po¬ 
larization, and position angle of the linear polarization. The triangles are 
points with no errors estimation. 


spectrograph was set to operate with 600 1/mm VPH grating, 1.68 
arcsec slit and the GG 385 blocking filter, yielding a spectral res¬ 
olution of 7 A FWHM in the range from 4350 to 7005 A. Three 
1200 s exposures were obtained, reaching individual S/N ratio of 
45 at the continuum. Quartz lamps calibration flats and bias im¬ 
ages were taken to correct the CCD signature, while CuHeAr lamp 
exposures were used for wavelength calibration, which resulted in 
typical 0.8 A or about 45 km s _1 calibration RMS residuals. The 
[OI] 5577 A telluric spectral line was used to assess the calibration 
accuracy. Two spectra of the HR 1544 spectrophotometric standard 
dHamuvIl 19921) were used for flux calibration. The data reduction, 
spectra extraction and calibrations were performed with IRAF stan¬ 
dard routines. The average spectrum of MLS110213 is presented 
in Fig. [3] 
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Wavelength (A) 

Figure 3. Average spectrum of MLS110213. 


3 RESULTS 

3.1 Light curves and ephemeris 

We analysed the photometric Ic, Rc and white light data obtained 
at OPD as well as the V data available from CRTS. This later 
dataset consists of 513 data points from the CSS and MLS sur¬ 
veys, spanning 8.3 years with a typical cadence of 4 consecutive 
exposures every visit, which are typically separated by 2-30 days. 
CRTS data points have magnitudes in the range from 16.5 to 19.5 
mag and are shown in Fig. [4] 

MLS110213 light curves obtained at OPD have higher tempo¬ 
ral resolution than CRTS data and clearly show eclipses, whose am¬ 
plitude reaches ~3 mag. We determined the timings of the six OPD 
eclipses as the midpoints of the bottom-of-eclipse light curves, the 
uncertainty being half of the exposure plus readout times (TableQ. 
Our measurements were converted to BJD(TDB) using the online 
code developed by lEastman. Siverd. & Gaudil d2010l) . These tim¬ 
ings were used to a first estimate of MLS110213 ephemeris. We 
then used the CRTS data to visually refine the ephemeris. The re¬ 
sulting ephemeris is 

T ecl = (BJD )2 455 862.6002 (±3) + 0.15779878(±6) xE.( 1) 

Errors of To and of the period are given in brackets. Fig. [5] 
shows the OPD white light, Ic and Rc light curves as well as 
CRTS data, folded with the period and epoch given in the above 
ephemeris. Despite the poor temporal resolution and the noise, 
long-spanning CRTS data are very well described by the ephemeris, 
with the eclipses clearly visible both in the high and low states (see 
below) with correct phasing. The eclipse has a flat bottom, with a 
duration of about 15 min. Its ingress is steep and lasts less than 2 
min, while the egress is less sharp. The shape of the eclipse is con¬ 
sistent with the eclipsed object being the WD and/or a post-shock 
region. 

By inspecting the CRTS light curve we note that it can be di¬ 
vided in two distinct brightness levels (high and low) separated by 
abo ut 1.5 mag, as can be seen in Fig. [4] (top and middle). Follow- 
in g IWu & Kiss! l l2008lt . we used a histogram of the magnitudes to 
choose the limit magnitude between high and low states. Figure [4] 
(bottom) shows a clear minimum at 17.85 mag, which we adopted 
as this limit. However, since MLS 110213 is an eclipsing system, 
the classification of a point as low or high can be tricky. We per¬ 
formed the classification of datasets of consecutive points (a visit) 


Table 2. Mid-eclipse timings of MLS 110213. 


E 

(cycles) 

T min 

(BJD(TDB) 2450000+) 

-468 

5788.750530(±72) 

-221 

5827.726780(± 129) 

0 

5862.600061 (± 118) 

7 

5863.704851 (±73) 

1845 

6153.738822(±75) 

4222 

6528.826641 (d= 16) 


as described below. The vast majority of the visits spans more than 
15 min, which is the eclipses length. Therefore most visits have at 
least one point out of the eclipse, and that point was used to esti¬ 
mate the brightness state of the corresponding dataset: if at least one 
point of a visit was brighter than 17.85 mag, then we considered this 
dataset as belonging to the high state. On the other hand, datasets 
with out-of-eclipse points fainter than 17.85 mag were classified as 
in the low state. But some CRTS visits are difficult to classify as 
high or low state because they include out-of-eclipse points consis¬ 
tent with both states. We tried to use a different limiting magnitude, 
but it did not help. Some of these points are seen around orbital 
phases 0.3 and 0.8 in Fig.[4](middle). A possible explanation is the 
presence of clouds during CRTS observations. 

Fig. 0 (top) shows MLS 110213 light curve in HJD. There is 
no typical duration of the high or low states, being the transition be¬ 
tween states fast, lasting about a day, which is a common behaviour 
among polars. According to this classification in distinct brightness 
levels, all of our OPD photometric data seem to correspond to the 
high state, except for I c band from August 14, 2011. 

The out-of-eclipse light curves show orbital variation in all 
bands, but it seems to depend on the brightness state (Figs. Q] and 
0. The CRTS data in low state have a modulation of around 1 mag 
with a minimum around 0.5 orbital phase. The phase dispersion is 
around 0.5 mag. In high state, the light curve has a larger disper¬ 
sion, of about 1 mag, and no orbital modulation is observed. More¬ 
over, possible departure from the high state may also be present 
in OPD data, as can be seen in the Ic light curves in Fig. |T] which 
were obtained one year apart. We did a careful analysis of objects in 
the field, and we discarded variation in the comparison star. Flick¬ 
ering with average 0.02 mag amplitude and time-scales of few min¬ 
utes is also visible in the out-of-eclipse OPD light curves. 

3.2 Polarization 

The polarimetric light curve (Fig. [2} confirms MLS110213 classifi¬ 
cation as a polar, since it shows variable polarized emission. These 
data were probably obtained while it was in the high brightness 
state because the CRTS data collected near our polarimetric obser¬ 
vations are also in high state. 

The system presents circular polarization varying in the range 
— 13 to 4 per cent in Rc band. The circular polarization does not 
have a standstill at zero, which indicates that a cyclotron compo¬ 
nent is always contributing to the observed flux in the high state. 
There are inversions of the polarization sign around phases 0.5 and 
0.8. This phenomenon indicates either one emitting region always 
visible, but observed from the backside in some phases, or two ac¬ 
cretion regions. The polarization has its maximum absolute value 
around phase 0.0, which corresponds approximately to the phase of 
maximum flux in the out-of eclipse low state light curve. It could 
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Figure 4. Top panel: V band light curves of MLS 110213 showing the 
high (open circles) and low (filled circles) brightness states along 8.3 years 
of CRTS observations. Middle panel: The same data as in the top panel, 
folded using our ephemeris. Bottom panel: Histogram of the magnitudes of 
MLS110213 using a bin size of 0.1 mag. 

be interpreted as if the low state modulation is also mainly due to 
cyclotron emission from the accretion column. Moreover, both in 
high and low states, the flux modulation is neither consistent with 
ellipsoidal variation nor with a heated secondary surface. 

The linear polarization is small, below 5 per cent, and vari¬ 
able. Its maximum values occurs at phases 0.4 and 0.75, which 
are near the phases when circular polarization changes sign: this is 
consistent with the direction of the emitting region changing from 
forward to backward. However, our data do not have a good SNR, 
and these conclusions should be taken with caution. 



Orbital phase 


Figure 5. Light curves of MLS 110213 folded with the ephemeris estimated 
in Section im Light curves expanded in the phase interval 0.9-1.1 to detail 
the eclipse shape. From top to bottom: OPD Re band, OPD Ic band, OPD 
white light, CRTS high state V band and CRTS low state V band. White 
light magnitudes are not calibrated. The blue dot-dashed line in the white 
light data (middle panel) is the result of the eclipse fitting discussed in Sec¬ 
tion^ 


3.3 Spectroscopy 

The average spectrum of MLS110213 is typical of polars (Fig. [3}. 
It is dominated by the Balmer and He 11 4686 A emission features, 
and also presents several lines of He I and the Bowen CIII/N III 
complex at 4630-4660 A. Emission lines of Fell 5169 A and Sill 
6347, 6371 A are also visible. The continuum is flat with a very 
slight blue slope and the Balmer lines have inverted decrement 
which, together with the intense He II 4686 A line, are characteris¬ 
tic of high-ionized optically-thick emission regions. Identification, 
equivalent width and FWHM of the main emission features in the 
average spectrum are presented in Table[3] 

The spectrum shows no sign of stellar features. Compar¬ 
ing the flux level to the CRTS data, we see that our spectra of 
MLS110213 were probably obtained in high brightness state. In 
addition, the absence of cyclotron humps and of secondary features 
in the spectra is consistent with the system being in a high accretion 
state, since these features may be diluted by the strong accretion ra¬ 
diation. 

Our spectra do not cover a complete orbital cycle. In the anal¬ 
ysis of each of the three individual spectra (Fig. 0, the He II 4686 
A and Balmer lines have profiles changing from symmetric, at the 
first spectrum, to profiles with extended red wings in the third spec¬ 
trum, taken 40 min later. We calculated the radial velocities of 
the C III/N III complex, He II 4686 A and Balmer lines, as a func¬ 
tion of the orbital ephemeris, by fitting a gaussian function to the 
peak of the lines (Table |4j- The lower limit to the RV amplitude is 
300 km s -1 . As these velocities were measured at the line peaks, 
they should map the velocities of the narrow component seen in 
polars, probably associated with the mass flux close to the donor 
star. 
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Figure 6. Individual spectra of MLS 110213 obtained at orbital phases 0.66, 
0.75 and 0.84 (from bottom to top, respectively). The spectra were shifted 
in the vertical axis for clarity, maintaining relative flux ratios. 


Table 3. Identification. EW and FWHM of emission features in the average 
spectrum of MLS110213. The FWHM measurements were obtained from 
Voigt profile fitting. 


Identif. 

—EW (A) 

FWHM (km s -1 ) 

He I4388 A 

4 

750 

He 14471 A 

16 

1070 

Clll/N 111 4640 A 

14 

1100 

He II 4686 A 

49 

640 

Hfl 

70 

680 

He I 4922 A 

9 

550 

He 15016 A 

13 

600 

He II 5412 A 

10 

610 

He I 5876 A 

21 

560 

Ho 

80 

590 

He I 6678 A 

12 

540 


4 CYCLOTRON MODELS 

We performed a modelling of the cyclotron emission of 
MLS110213 . CYCLOPS code was used to fit the available photo¬ 
metric and polarimetric data and to estimate some geometrical and 
physical properties of the system. CYCLOPS is a radiative trans¬ 
fer code that calculates the cyclotron and bremsstrahlung emission 
from a 3D post-shock region. The accretion geometry is defined 
by a bipolar magnetic field, B po i e , centred in the WD. The plasma 
temperature ( T ) and density (N e ) structures are defined by shock¬ 
like profiles in the radial direction and by an exponential decay or 

Table 4. Radial velocities of selected lines, in km s 1 . 


phase 

C iii/N III 4640 A 

He il 4686 A 

h /3 

H a 

0.66 

-124 

-106 

-129 

-127 

0.75 

25 

19 

-12 

-30 

0.84 

186 

169 

156 

98 


a constant profile in the tangential direction. The shock-like pro¬ 
files are analytical function s having a sha pe s imilar to the magne - 
tised ID shock solution of iCroppeii 1199ol) an d lSaxton et al.1 1200 7). 
Their parameters are the maximum temperature, the maximum den¬ 
sity, and the height of the post-shock region. These functions are 
presented in equations (1) and (2) of lSilvaet alj < 120131) . 

Eleven parameters define a model, as listed in the upper part of 
Table[5] The parameters define the geometrical and physical prop¬ 
erties of the the post-shock region, which is divided into voxels. 
The spatial resolution of a model is a free parameter, which de¬ 
fines the voxel size. The emissivities of the four Stokes parameters 
are calculated according to the B, T, and N e of each voxel, and 
after that, the radiative transfer is calculated, from the bottom to 
the top of the region, according with the line of sight. The result 
is an image in each orbital phase for each Stokes parameter. These 
values are integrated in area to obtain the flux and polarization as a 
function of the orbital phase. The figure-of-merit to quantify the fit¬ 
ting quality is the ,\ 2 . The genetic algorithm PIKAIA l lCharbonneaul 
119951) is used to identify the regions in the parameter space having 
the best models, and then an amoeba based code refines the search. 
The model can fit several optical bands simultaneousl y . For a de- 
tailed desc r iption of the code, see ICo sta & Rodrigues) d2009h and 
ISilva et alJ d2013h . 

We fitted the high state data of MLS110213 comprising V 
light curve from CRTS, R c light and polarization curves from OPD 
and I c light curve obtained 2012 August 13, also from OPD (Fig. [7]). 
Magnitudes and polarizations were converted into total and polar¬ 
ized fluxes. The data in each band were binned in 40 phase inter¬ 
vals. Although the data are not simultaneous, we consider that the 
brightness states are similar. The large error bars of the V light 
curve are probably the result of a combination of measurements 
taken in slightly different brightness states. Photometric measure¬ 
ments obtained near eclipse, between phases 0.95 and 0.1, were 
analysed separately (see Section O, as CYCLOPS does not handle 
the eclipse of the post-shock region by the secondary. 

In our minimisation procedure, we keep free all parameters 
listed in Table 5 except for the inclination, which was limited be¬ 
tween 75 and 90 deg due to the presence of total eclipse. In the 
modelling, we used the following frequencies to represent each op¬ 
tical band: V = 5.45 xlO 14 Hz, Rc = 4.49xlO 14 Hz and Ic = 
3.8xl0 14 Hz. The results of our modelling are presented in the 
next section. 

4.1 Modelling discussion 

As already discussed in Section [3~2l the inversion of the polariza¬ 
tion sign may be observed in two situations: (i) two emitting re¬ 
gions; (ii) one emitting region whose magnetic field lines alternate 
direction back and forth along the orbital cycle. As we found a rea¬ 
sonable fit using only one region, we adopted this solution, there¬ 
fore minimising the number of free parameters of the model. 

Using PIKAIA with distinct parameters intervals, the best fits 
were always found in a same group of models. The fits considering 
layers with constant values of N e and T are worse than those in 
which these quantities are let to vary in the direction perpendicular 
to the column height. The model having the smallest \ 2 is shown 
in Fig. [7] and its parameters are presented in Table [5] We discuss 
this model below. 

The model reproduces quite well the Rc circular polarization 
(Fig. [7] left/top panel, black lines). In particular, the observed sign 
inversion is well fit. The observed linear polarization flux is small 
and noisy, but larger than that produced by the model. 
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Figure 7. MLS 110213 data (points and triangles- the last with no erros estimates) and the best-tit model (lines) described in text and whose parameters are 
presented in Table[5] Left panel: Polarization modelling. From top to bottom, the panels show circular polarized flux ( Fy ), linear polarized flux (Fp) and 
angle of linear polarization ( 6 ). The black dots are the Rq polarization data and the black solid line represents the best Tic model. The blue dotted line is the 
V band model prediction and the red dashed line represents the model prediction in Ic band. Right panel: Total flux modelling (Ft). From top to bottom, the 
panels shows the light curves in V, Rc and Ic bands. The black lines are the best models. 


The flux modulation in V, Rc, and Ic bands are reasonable 
well fit (Fig. [7] right panel). Figure [7]also presents the polarization 
predictions to V and Rc bands (blue dotted and red dashed line, 
respectively). The predicted circular polarization in V band is sim¬ 
ilar in shape, but has a larger amplitude than Rc band, and the sign 
inversion is more pronounced. I c band model prediction shows null 
circular polarization from phase 0.5 to 1.1. 

Figure[8]is a geometrical cartoon of our modelling of the emit¬ 
ting region on the WD surface along the orbital cycle. From phase 
0.7-0.1 the emitting region is visualised from the top, when the cir¬ 
cular polarization is negative and larger. From 0.2 to 0.6 the region 
is partially occulted and is seen from the opposite direction, hence 
the circular polarization is positive and smaller. The sign inversions 
can be associated with the following configurations: the first sign 
inversion occurs around phase 0.5, when the region is seen com¬ 
pletely by the opposite side for the first time; the second inversion 
occurs around phase 0.8, when the entire region starts to point to 
the observer again. Therefore, we can conclude that the circular 
polarization sign inversion can be due the changes in the angle of 
view of one dominant emitting region. The change in the signal 
does not happen early in phases 0.2 or 0.3 because of the partial 
self-eclipse. The self-eclipse could also decrease the value of the 
absolute circular polarization flux. 

The best fit model corresponds to an extended region, with a 
relatively high magnetic latitude and moderate electronic tempera¬ 
ture. Orbital inclination from 75° to 81° produces very similar \ 2 
values: inclinations out of the range 74° - 85° increase \ 2 by at 
least 30%. The best-fit magnetic field in the magnetic pole, B po i e , 
is 35 MG. This corresponds to a magnetic field in the region varying 
from 18 to 33 MG. Values of B po i e between 34 and 36 MG produce 
very similar \ 2 values, while B po i e of 30 and 40 MG produce \ 2 
30% larger. The region mean temperature is 9 keV. The mean tem¬ 


perature weighted with the square density, T pon d, is around 5 keV. 
For the temperature and density we found that the ranges 24 - 30 
K and 14 - 14.3 cm^ 3 provide small changes in the \ 2 - Out of the 
limits of 22 - 35 K and 13.8 - 14.6 cm~ 3 , the \' 2 value increases 
by at least 30%. These are low values of temperature for a polar 
post-shock region and indicate a dominance of soft emission in X- 
ray. The maximum temperature in the emitting region, T max, corre- 
sponds to the temp e rature in the sh ock front. U sing the approach of 
iFuiimoto & Ishidal d 19971) and the iNauenbergl d 19721) mass-radius 
relation, we estimated a mass of 0.67 Mg for the WD from T max- 
The orbital period of MLS1T0213 translate s to a secondary mass, 
M 2 , of 0.279 Mg dKnigge. Baraffe & PattersoiJl201 ll) . resulting in 
a mass ratio, q, of 0.42. The stellar parameters are also summarised 
in Table 5. 

The presented modelling of the accretion column can be re¬ 
fined using new additional data, for instance, observations in other 
broad bands, simultaneous data acquisition and/or spectropolarime- 
try is desirable. 


5 ECLIPSE FITTING 

We fitted the eclipse profile seen in the clear band observations 
considering a simple model of two spherical stars in a circular orbit. 
The clear data present the best temporal resolution covering the 
eclipse. The eclipse is box shaped indicating a total eclipse of a 
very small object. 

We determined M2=0.279 Mg and R. 2 = 0.37 Rg by 
interpolating the semi-empirical CV donor star sequence 
dKnigge. Baraffe & Patterson! 1201 ll) as a function of the or¬ 
bital period of the system. The flux of the donor star is considered 
negligible and Mi is obtained from the mass ratio, q = M 2 /M 1 . We 
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Table 5. Physical and geometrical properties estimated for MLS 110213. 


Cyclops input parameters 

Fitted values 

i 

n° 

p 

39° 

A long 

87° 

A.r 

0.29 

h 

0.18 Rwd 

Si 

0.29 

B p ole 

35 MG 

Blat 

65° 

Blong 

490 

Tmax 

29 keV 

Ne max 

14.3 cm - 3 

Model results 

Values 

Breg 

18-33 MG 

(T) 

9 keV 

r p 

- L pona 

5 keV 

Trange 

1-29 keV 

Sphase 

- 0.11 

x 2 

0.395 

Stellar parameters 

Values 

<7 

0.42 

Mi 

0.67 Mq 

Ri 

0.012 R 0 

m 2 

0.28 M 0 

r 2 

0.36 R 0 


Parameters description - i : orbital inclination; (3: colatitude of the reference 
point of the coordinate system inside the emitting region; A i ong : one-half 
of the azimuth of the threading region; A#: one-half of the radial exten¬ 
sion of the threading region; h : height of the post-shock region; fi : position 
in the longitudinal direction of the threading point relative to the geomet¬ 
rical centre of the threading region; B po i e : polar magnetic field intensity; 
Bi a f : latitude of the magnetic axis; Bi ong \ longitude of the magnetic axis; 
Tmax ■ maximum electronic temperature; N e max : log of maximum elec¬ 
tronic density; B reg \ magnetic field in the post-shock region; (T): mean 
eletronic temperature; T pon( i: mean temperature weighted by the square 
density (^( T vox N% ox )/ ^ N^ ox ;T ra nge'- range of temperatures of in 
the post-shock region; S p h as e'■ phase shift applied to the model. 


use the masses to determine th e dist ance between the stars from 
the Kepler Third Law and use iNauenberel dl972h relationship to 
determine Ri from a given Mi. Thus the parameters of the eclipse 
fitting are the orbital inclination i and the mass ratio q. 

There is a degeneracy in the possible values of i and q that can 
fit the eclipse profile. Mi will reach the limit of 1.4 Mg at q = 0.20 
and i = 19°. We can increase q while decreasing i until it reaches 
the minimum inclination limit for a total eclipse at i = 73° and q 
= 0.75. Therefore, the observed eclipse width can be reproduced 
using the following intervals of parameters: Mi = 0.37 - 1.4 Mg, 
with q = 0.75 - 0.20 and i = 73° - 19°. 

The CYCLOPS modelling shows no significant variation in the 
fitting for the interval i = 15° - 81°, indicating that the CYCLOPS 
geometry can also describe the observed eclipse profile. The blue 
dot-dashed line in the middle panel of Figure [5] shows the eclipse 
fitting using q = 0.42, as estimated in the previous section. This 
value corresponds to Mi= 0.67 Mg and Ri = 0.012 R 0 , keeping 
M 2 constant. The inclination that fits the eclipse profile is i = 15.8°, 
in agreement with CYCLOPS estimated interval. 

In polars, the emitting region can be much brighter than 
the WD itself, therefore determining the eclipse profile. The 
MLS 110213 eclipsed region is much smaller than the secondary 
star, as shown by the fast ingress and egress of the total eclipse. But 
the data do not have enough time resolution to allow the determi¬ 
nation of the size of the eclipsed region. 

We consider that the approximation of a spherical sec¬ 
ondary star is valid because, at conjunction, the projected 
area of the secondary star on the plane of sky is min¬ 
imum. To check this assumption, we com pare our results 
with those of Ichanan, Middleditch. & Nelsonl < il976l) . who con¬ 
sider the actual geometry of a Roche-lobe-filling companion 
and a pointlike eclipsed object. We measured the half-eclipse 
width of MLS 110213 as 9.36°. Accordingly with Table 1 of 
IChanan, Middleditch, & Nelsonl d 1976h . the mass ratio of q=0.42 
results in an i nclination of i = 16.1° , in ag reement with our re¬ 
sult. Moreover. 1 Andronov & Andrvchl d2014l ) compare the use of a 
spherical approximation and an elliptical projection onto the celes¬ 
tial sphere for a secondary that fills its Roche lobe and found that, 
in the spherical approximation, the inclination is underestimated by 
about 2°, a value which is inside our uncertainties. 




0.5 0.6 0.7 0.8 0.9 

Orbital Phase 


Figure 8. Geometrical representation of the emitting region on the WD sur¬ 
face in MLS 110213. The post-shock region is represented by its walls (blue 
lines). The curved red line near the centre of the emitting region is a mag¬ 
netic field line threading the emitting region. It is shown to represent the 
accretion column geometry. The radial green line is the magnetic axis. 


6 DISTANCE ESTIMATION 

In order to estimate the distance to MLS110213, we adopted 
the Period-Luminosity-Colo urs (PLCs) method derived by 
lOzddnmez. Ak & Bilin 1 120 151) The method uses the orbital period 
of the s ystem together with J , K, and W 1 magnitudes, from 
2MASS dSkrutskie et al.l 1 20061) and WISE J Wright et al.l l20ldl . 
to estimate abso lute magn itudes of CVs. The method is based on 
the iBames & Evansl dl976h relation which states that, for late-type 
main sequence stars, the surface brightness is nearly constant in the 
near-infrared. The PLCs relation is calibrated by the trigonometric 
parallaxes of 25 CVs. The accuracy in Mj provided by the PLCs 
method is ±0.29 mag. As the PLCs method does not consider 
other sources of infrared light except the secondary star, it yields a 
lower limit to the distance. The PLCs relation is: 

Mj = 5.966 - 4.781 log P + 5.037( J-K s )o + 0M7(K S -Wl) 0 (2) 

where P is the orbital period given in hours and the de-reddened 
colours (J — Ks) 0 and ( Ks — WT)o are calculated from 
Jo = J — Aj, Kso = Ks — Ak s and Wlo = Wl — Awi - We 
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adopted Aj = 0.887 x E(B - V), A Ke = 0. 382 x E(B - V) 
and Aw 1 = 0.158 x E(B — V) dOzdonmez, Ak & Bilin 20151) . 

The infrared magnitudes of MLS 110213, obtained from the 
2MASS and WISE database, are J = 15.928, H = 15.264, 
K s = 15.000, I'Ll = 14.680 and W2 = 14.163. The 

colour excess in the direct ion of MLS 110213, obtained from 
ISchlaflv & Finkbeineri d201 ill maps, is E(B — V) = 0.0941. With 
those values and considering the accuracy in Mj, we estimate a dis¬ 
tance to MLS110213 of 406±54 pc, assuming that the secondary is 
the only infrared light source in the system. Despite being a lower 
limit, this distance should not be very different from the true value, 
since our CYCLOPS modelling indicates that the cyclotron contri¬ 
bution to the total infrared flux of the binary system is small. This 
may not be the case for other magnetic CVs. In non-magnetic CVs, 
the infrared flux may have contributions from other structures like 
the outer regions of the accretion disc. 


7 DISCUSSION AND CONCLUSIONS 

We have obtained photometric, polarimetric and spectroscopic data 
of MLS 110213, which is one of the polar candidates identified and 
monitored by the Catalina Real Time Transient Survey. We confirm 
the system as a new eclipsing polar. 

The optical spectrum, obtained in the 0.7 - 0.88 phase inter¬ 
val, is very similar to the spectra of other polars. It shows a flat 
continuum and an inverse Balmer decrement, and also conspicuous 
features of highly ionized optically thick regions, like strong emis¬ 
sion lines of He II 4686 A and CIII/N III 4640 A. There is no clear 
evidence of the secondary star in our spectrum. 

Photometric data obtained at CRTS and OPD show that 
MLS 110213 presents deep (3 mag) eclipses with flat bottoms, indi¬ 
cating that the system has high orbital inclination. From OPD mid¬ 
eclipse timings we estimated an orbital period of 3.787 h, which 
places MLS 110213 as the 5 th eclipsing polar found above the or¬ 
bital period gap. 

The long term observations obtained by CRTS clearly show 
the system in two distinct brightness states, separated by 1.5 mag. 
In the high brightness state the out-of-eclipse light curve is almost 
flat, while in low state a strong modulation is observed, with a max¬ 
imum at phase 0.175. The presence of distinct brightness states is 
typical of polars and are understood as variations at the mass accre¬ 
tion rate. 

MLS 110213 presents circular polarization in the R c band 
(obtained in high state) varying from -13 to 4 per cent with a max¬ 
imum around phases 0-0.125, consistent with cyclotron emission 
from the accretion column. The polarization does not have a stand¬ 
still at zero, and there are phases when the signal changes, which 
may indicate either one emitting region always visible, but ob¬ 
served from the backside in some phases, or two accretion regions. 

The CYCLOPS modelling of high state V, R c and I c light 
curves and of one R c polarimetric curve provided an estimate of 
the system parameters and a description of the main observed fea¬ 
tures, including the eclipse profile. The model agrees with the in¬ 
terpretation that MLS 110213 has one main circumpolar emitting 
region, although the possibility that it could be a two-pole accretor 
cannot yet be ruled out. This model constrains the magnetic field 
from 30 to 40 MG and the orbital inclination from 75° to 85°. The 
eclipse fitting corroborates CYCLOPS results. Nevertheless, addi¬ 
tional polarimetric data in V and I c band are necessary to a better 
determination of the system parameters. 

The field of MLS 110213 was observed by ROSAT, RXTE, 


Integral and Swift observatories, but the source was not detected 
by none of them. Using data from the JEM-X camera at Inte¬ 
gral, we estimated an upper limit of 1.3xlO~ 10 erg cnW 2 s _1 
for the X-ray flux in the energy range of 2-5 keV. This flux 
is quite high, and we consider that the Integral non-detection 
is due to a short integration time. The ROSAT Bright Survey 
detection limit flux is 2xl0 -12 erg cm -2 s -1 in the 0.12- 
2.48 keV band, consideri ng a ty pi cal m odel for X-ray polars emis¬ 
sion IPretorius. Knigge. & Schwopel2013h . Adopting this flux limit 
and the distance estimate of 406 pc, the X-ray luminosity of 
MLS110213 is lower than 3.9x 10“ 31 erg s _1 . All ROSAT polars 
with X-ray lumi nosities similar or smalle r than this value are nearer 
systems ( iPretorius. Knigge. & Schwopelboijl) . Therefore the low 
X-ray flux of MLS 110213 is possibly explained by the distance 
to the source and accentuated by the low temperature of the emit¬ 
ting region, as estimated by the CYCLOPS fit. This fit predicts a 
modulation of the X-ray light curves with the minimum flux from 
phases 0.2 to 0.5, caused by partial self-eclipse of the emitting re¬ 
gion. Deep X-ray observations of MLS 110213 would be very use¬ 
ful to constrain the geometric and physical properties of the system. 

MLS 110213 is among the few polars not found by X- 


ray surveys. Some 
rate polar (LARPs), 

of them are classified as low accretion- 
a class now also recognised as pre-polars 

(e.g., Schwo 

pe et al. 

2002|_ Webbink & Wickramasinghe 2005; 

Schmidt et al. 

2005; 

Breedt et alJ 20 ill). MLS 110213 shares the 


low X-ray flux and low temperature of this class, but the low X- 
ray brightness may be explained by the MLS 110213 distance. Its 
post-shock region seems to be dense, produced by a Roche lobe 
overflow mass transfer. This diverges from MLS 110213 being a 
pre-polar, since in them the mass transfer is supposed to originate 
from a secondary wind, with mass transfer rates about 100 times 
smaller than in standard polars. The most important signature of a 
low mass transfer rate, the cyclotron harmonics, is also lacking in 
theMLS110213 spectra. Therefore, we do not classify MLS 110213 
as a LARP/pre-polar. 

iLi & Wickramasinghel dl998i . and a series of previous arti¬ 
cles) study the effect of the WD magnetic field on the secondary 
magnetic braking. They show that the interaction of the magnetic 
fields of both stars can decrease the magnetic braking. This effect 
is more pronounced if the magnetic dipole axis is parallel to the 
WD spin axis. Hence they concluded that the evolution of mag¬ 
netic CVs should take into account the direction of magnetic dipole 
axis. They argued that most polars tend to have magnetic field dipo¬ 
lar axis approximately aligned with the WD spin axis (at angles 
smaller than 25°) and, in this situation, the MB can be lessened by 
at least an order oUmagnitude in compa rison with non-magnetic 
CVs. lWebbink & Wickramasinghel d2002i) shows that this decrease 
in the magnetic braking can populate the period gap in polars. The 
inclination of the WD dipolar magnetic field obtained by our CY¬ 
CLOPS modelling of MLS1 10213 is 25°. This value is consistent 
with the "aligned" field scenario. 

MLS 110213 is a new eclipsing polar above the period gap. Its 
post-shock region is cooler than in most mCVs, indicating a low- 
mass WD. MLS 110213 populates a rare sub-group of polars, near 
the upper limit of the period distribution, important to understand 
the evolution of mCVs. Our modelling is a first step to determine 
parameters as the magnetic field intensity and inclination, which 
could be used to constrain evolutionary models of mCVs. 
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